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Abstract 

The physics potential of an intense source of low-energy muons is studied. 
Such a source is a necessary stage towards building the neutrino factories and 
muon colliders which are being considered at present. The CERN Neutrino 
Factory could deliver muon beams with intensities 3-4 orders of magnitude 
higher than available now, with large freedom in the choice of the time struc- 
ture. Low-energy muon physics contributes to many fields of basic research, 
including rare muon decays, i.e., decays that do not conserve muon number, 
measurements of fundamental constants, the muon anomalous magnetic mo- 
ment, determination of the Lorentz structure of the weak interaction, QED 
tests, CPT tests, proton and nuclear charge distributions (even for short-lived 
isotopes), and condensed matter physics. In studying the experimental pro- 
gramme, we analyse the present limitations, list the requirements on the new 
muon beams, and describe some ideas on how to implement these beam lines 
in a CERN neutrino factory complex. 
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1 INTRODUCTION 



Ever since the discovery of the muon the study of its properties has contributed to a deeper understanding 
of Nature at the smallest distance scale. Muon physics played a fundamental role in establishing the V-A 
structure of weak interactions and the validity of quantum electrodynamics. Moreover, muon physics has 
not yet exhausted its potential and, indeed, may provide crucial information regarding one of the most 
fundamental quests in modern physics: the structure of the theory which lies beyond the Standard Model 
of particle physics. 

Muon storage rings are currently, and seriously, being considered as options for future endeav- 
ors of accelerator laboratories such as CERN and Fermilab. The primary aim of the machine presently 
considered is the study of neutrino properties, hence they are also referred to as neutrino factories (NU- 
FACT). They are a first step towards a still more ambitious goal: muon colliders, capable of 'precision' 
Higgs physics and the exploration of the highest energy frontier. 

Neutrino factories are also an ideal place to study muon properties, since they provide, necessarily, 
muon fluxes which are orders of magnitude larger than what can be obtained at present. It is, therefore, 
imperative to understand how to take full advantage of these intense muon beams in order to significantly 
improve on the reach of stopped-muon experiments. This is the aim of this part of the report. 

While precise measurements of the muon lifetime and Michel parameters provide tests for the 
theory of weak interactions and its possible extensions, one of the main interests in muon physics lies 
in the search for processes that violate muon number. The discovery of decays such as /i + — ► e + 7 and 
/i + — > e + e~e + or of \i~-e~ conversion in nuclei would be an indisputable proof of the existence of new 
dynamics beyond the Standard Model. 

Global symmetries (like individual lepton numbers), as opposed to local symmetries, are consid- 
ered not to be based on fundamental principles and are expected to be violated by gravitational effects, 
in the strong regime, and, more generally, by higher-dimensional effective operators which describe lo- 
cal interactions originating from some unknown high-energy dynamics. Baryon number conservation 
is another example of an Abelian global symmetry of the Standard Model, which can be broken by 
new-physics effects. 

Atmospheric and solar neutrino experiments have provided strong evidence for neutrino oscilla- 
tions. This implies violation of individual lepton numbers (Li) and, most likely, of total lepton number 
(L), which is a first indication of physics beyond the Standard Model. Current neutrino data indicate val- 
ues of the neutrino masses corresponding to non-renormalizable interactions at a scale M ~ io 9 ~ 14 GeV. 
New lepton-number violating dynamics at the scale M cannot yield observable rates for rare muon pro- 
cesses, since the corresponding effects are suppressed by (m^/M) 4 . The observation of muon-number 
violation in muon decays would thus require new physics beyond the one responsible for neutrino masses. 
Theoretically, however, there is no reason why Lj and L would be broken at the same energy scale. In- 
deed, in many frameworks, such as supersymmetry, the Lj breaking scale can be close to the weak scale. 
In this case, muon processes with violation would occur with rates close to the current experimental 
bounds. 

It is also very important to stress that the information which can be extracted from the study of 
rare muon processes is, in many cases, not accessible to high-energy colliders. Take supersymmetry as 
an example. While the LHC can significantly probe slepton masses, it cannot compete with muon-decay 
experiments in constraining the slepton mixing angles. 

In Section || we discuss 'normal' muon decay and its importance to fundamental physics. In 
Section 3 we address muon-number violating processes, first discussing the theoretical expectations and 
then addressing the issue of how to improve on current experimental bounds given a larger muon flux. 
In particular we discuss which elements of the experimental set-ups need to be improved or drastically 
modified. In Section 4 we address the measurements of fundamental properties of the muon, and in 
Section 5 we discuss bound muon systems (muonic atoms and condensed matter systems). In Section 6 
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we summarize the beam requirements imposed by the various experiments and discuss how these could 
be met in a neutrino factory complex. Our conclusions are given in Section 7. 

2 NORMAL MUON DECAY 

2.1 Theoretical background 

All measurements in normal muon decay, n~ — ► e~V e v^, and its inverse, z^e - — ► fi~u e , are successfully 
described by the ' V-A' interaction, which is a particular case of the local, derivative-free, lepton-number- 
conserving, four-fermion interaction The 'V-A form and the nature of the neutrinos (z7 e and u e ) have 
been determined by experiment [Q, ||]. 

The observables in muon decay (energy spectra, polarizations and angular distributions) and in 
inverse muon decay (the reaction cross section) at energies well below m^yc 2 may be parametrized in 
terms of the dimensionless coupling constants g2p and the Fermi coupling constant Gf. The matrix 
element is 

M = ^ g2»{ee\^\(v e )n){(v,)m\T,\^) ■ (D 

7 =S,V,T 

We use here the notation of Fetscher et al. , [|2|, ^ who in turn use the sign conventions and definitions of 
Scheck Q|. Here 7 = S, V, T indicate a (Lorentz) scalar, vector, or tensor interaction, and the chirality 
of the electron or muon (right- or left-handed) is labelled by e, /i = R, L. The chiralities n and m 
of the v e and the are determined by given values of 7, e and fi. The 10 complex amplitudes g]^ 
and Gf constitute 19 independent parameters to be determined by experiment. The 'V-A interaction 
corresponds to g^ L = 1, with all other amplitudes being 0. 

With the deduction from experiments that the interaction is predominantly of the vector type and 
left-handed [g^ L > 0.96 (90 %CL)], there remain several separate routes of investigation of normal 
muon decay, which will be discussed in the following. 

2.2 Muon lifetime measurements 

The measurement of the muon lifetime yields the most precise determination of the Fermi coupling 
constant Gf, which is presently known with a relative precision of 9 x 10~ 6 . Improving this measurement 
is certainly an interesting goal [gj, since Gp is one of the fundamental parameters of the Standard Model. 
A clean beam pulse structure with very good suppression of particles between pulses is indispensable. 
Presently three experiments are in progress, two of which are at PSI [Q] and one is located at RAL [||]. 
An improvement in the precision of 7^ by about a factor of 20 can be expected. An additional order 
of magnitude could be gained at a neutrino factory primarily from increased muon flux with the major 
systematics being pile-up and detector timing stability. 

There are two caveats, however: 

• Reducing the error on Gf by precise measurements of the muon lifetime would unfortunately not 
improve the electroweak fits, because the error on the dimensionless input GfM§ is dominated by 
the uncertainty on M§, which is now 48 x 10" 6 [M z = (91 188.2 ± 2.2] MeV [§]). 

• Gf is commonly determined assuming exclusively V-A interactions. A somewhat more general 
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formula has been given by Greub et al. [10]: 
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Here, besides the muon lifetime and the muon mass, radiative corrections to first order and mass 
terms are included. Most important is the muon decay parameter r\ which is in the SM. If 
we assume that only one additional interaction contributes to muon decay, then 77 ~ ^Reg RR , 
where g RR corresponds to a scalar coupling with right-handed charged leptons. Including the 
experimental value of 77 = (— 7 ± 13) x 10 -3 [ |TT| ] the error on Gp increases by a factor of 20. 



2.3 Precision measurement of the Michel parameters 

The measurement of individual decay parameters alone generally does not give conclusive information 
about the decay interaction owing to the many different couplings and interference terms. An example 
is the spectrum Michel parameter g. A precise measurement yielding the V-A value of 3/4 by no means 
establishes the V-A interaction. In fact, any interaction consisting of an arbitrary combination of g^ L , 
9hR' 9rL' 9rR' 9rr an( ^ 9ll wm yi e ld exactly g = 3/4 [|l2|]. This can be seen if we write g in the form 

m 



3 3 fi V i2 1 1 V |2l . r, (\ T i2 1 1 T i2\ , D / S T* . S T* \ ,c\ 

6- 4 = ~Z \\9lr\ +\9rl\ } + 2(|5lrI +|5rlI J + Re (<7lr<7lr + 9kl9kl ) ■ (5) 



For g = I and g RL = g^R = ( no tensor interaction) one finds g RL = g^R = 0> with au °f the 
remaining six couplings being arbitrary! On the other hand, any deviation from the canonical value 
certainly would signify new physics. Tree-level new physics contributions to the Michel parameters 
occur in supersymmetric theories with i?-parity violation or theories with left-right symmetric gauge 
groups. For instance, the i?-parity violating interactions AsnL^ Ll -^R + ^^^h^h^R^ (where 



the index denotes the lepton generation) give the following contributions [14] 
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For a left-right model, one finds 
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where flyy^ is the Wl-Wr mixing angle, and M\y 1 (Mw 2 ) is the mass of the mainly left (right) charged 
gauge boson. Measurements of g and £ with a precision of 10~ 4 can probe Wr masses of about 1 TeV (in 
the most unfavourable case $w R = 0) and values of the i?-parity violating couplings A311 ~ A322 ~ 0.2 
(for a slepton mass of 200 GeV). These tests are competitive with direct searches at high-energy colliders. 

There exist also observables which yield valuable information even if they assume their canonical 
values, all of which are related to the spin variables of the muon and the electron: 
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A measurement of the decay asymmetry yields the parameters 5 and P M £. Especially interesting 
is the combination P^S/g, which has been measured at TRIUMF [15] with a precision of « 



3 x 10 . A new, ambitious experiment measuring g, 5 and P„£ is currently under construction at 



TRIUMF [16] 



A measurement of the longitudinal polarisation of the decay electrons Pl consistent with 1 yields 
limits for all five couplings where the electrons are right-handed. This is a difficult experiment 
due to the lack of highly polarised electron targets used as analysers. The present precision is 
AP L = 45 x 1(T 3 . 

The angular dependence of the longitudinal polarisation of decay positrons at the endpoint energy 



is currently being measured at PSI by the Louvain-la-Neuve-PSI-ETH Zurich Collaboration Q17| ] 



This yields the parameter £" which is sensitive to the right-handed vector and the tensor currents. 

A measurement of the transverse polarization of the decay positrons requires a highly polarised 
pulsed muon beam. From the energy dependence of the component Pj^ 1 one can deduce the low- 
energy decay parameter rj which is needed for a model-independent value of the Fermi coupling 
constant. The second component Pr 2 , which is transverse to the positron momentum and the muon 
polarisation, is non-invariant under time reversal. A second generation experiment performed at 
PSI by the ETH Zurich-Cracow-PSI Collaboration has just finished data-taking [|TJ, |i"9[ ]. The 
expected experimental errors are APr^ = APt 2 = 5 x 10~ 3 . 



2.4 Experimental prospects 

As mentioned above, the precision on the muon lifetime can presumably be increased over the ongoing 
measurements by one order of magnitude. Improvement in measurements of the decay parameters seems 
more difficult. Most ambitious is the TRIUMF project which will soon start taking data and thereby 
improve the parameters g (positron energy spectrum), P M £ and 5 (decay asymmetry) by more than one 
order of magnitude. The limits on most other observables are not given by the muon rates which usually 
are high enough already («3x 10 8 s" 1 at the //El beam at PSI, for example), but rather by effects like 
positron depolarisation in matter or by the small available polarisation (< 7%) of the electron targets 
used as analysers. The measurement of the transverse positron polarisation might be improved with a 
smaller phase space (lateral beam dimension of a few millimetres or better). This experiment needs a 
pulsed beam with high polarisation. 



3 SEARCH FOR MUON NUMBER VIOLATION 



3.1 Theoretical considerations 

In the Standard Model (SM), muon number is exactly conserved. When neutrino masses are added and 
neutrino oscillations take place, muon-number violating processes involving charged leptons become 
possible as well. However, because of the smallness of neutrino masses, the rates for these processes are 
unobservable p^ , p| ] ; for instance 
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The observation of muon number violation in charged muon decay would, therefore, serve as an unam- 
biguous sign of new physics and indeed, a number of SM extensions may be probed sensitively by the 
study of rare muon decays. Here we will concentrate on supersymmetric models and models with extra 
dimensions, but it should be pointed out that various other SM extensions also predict observable rates 



for the rare p decays, like models with new Z' gauge bosons [22], leptoquarks [ ]23| ] or Lorentz-invariance 
violation fl24|, [25|]. For a recent review on muon number violation, see Ref. [53]. 
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We first present a model-independent formalism for studying rare muon decays and for compar- 
ing the rates of the different muon-number-violating channels. Then, we discuss predictions for the 
branching ratios of rare muon processes in extensions of the SM with low-energy supersymmetry, in the 
context of the seesaw mechanism and grand unification, both with and without the conservation of R 
parity. Finally, we present the expectations in models with extra dimensions. 



3.1.1 Model-independent analysis of rare muon processes 

Although a purely model-independent analysis based on effective operators cannot make any prediction 
for the absolute rate of rare muon processes, it can be very useful in determining the relative rates. We 
will compare the rates for fi + — ► e + 7, fi + — ► e + e~e + , and pT-e~ conversion. In a large class of 
models, the dominant source of individual lepton number violation comes from a flavour non-diagonal 
magnetic-moment transition. Let us therefore consider the effective operator 



+ h.c. 



(9) 



This interaction leads to the following results for the branching ratios of /i + — ► e + 7 (B(p — ► ej)) and 
/i + — ► e + e~e + (B(fi — > 3e)), and for the rate of pT-e~ conversion in nuclei normalised to the nuclear 
capture rate (B(/jN — ► eN))\ 
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Here B(A, Z) is an effective nuclear coefficient which is of order 1 for elements heavier than aluminium 
[p7|]. The logarithm in Eq. (11) is an enhancement factor for B(/j, — > 3e), which is a consequence of 
the collinear divergence of the electron-positron pair in the m e — ► limit. Nevertheless, because of the 
smaller phase space and extra power of a, B{p, — > 3e) and B(fiN — ► eiV) turn out to be suppressed with 
respect to — > &f) by factors of 6 x 10~ 3 and (2-4) x 10" 3 , respectively. This should be compared 
with the experimental sensitivities for the different processes, discussed in the next subsections. 

Next, let us include an effective four-fermion operator which violates individual lepton number 



£ = -jjiiin^zLhiiifL + h.c. 



(13) 



where / is a generic quark or lepton. The choice of the operator in Eq. ( fl3[ ) is made for concreteness, and 
our results do not depend significantly on the specific chiral structure of the operator. First we consider 
the case in which / is neither an electron nor a light quark, and therefore /i + — ► e + e~e + and \i~-e~ 
conversion occur only at the loop level. Comparing the p, + — ► e + 7 rate in Eq. (|10|) with the contributions 
from the four-fermion operator to B([i — > 3e) and B(fiN — ► eN), we find 
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Here is the number of colours of the fermion / and M f is the heavy-paiticle mass generating the 
effective operators (typically Mp is much smaller than A or Aj? because of loop factors and mixing 
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angles). The logarithms in Eqs. (14) and (15) correspond to the anomalous dimension mixing of the 
operator in Eq. (|l3|) with the four-fermion operator generating the relevant rare muon process [28]. If 
A ~ Ap, then the contributions from the four-fermion operator are irrelevant, since the ratios in Eqs. (11) 
and ([l2|) are larger than those in Eqs. ( |l4| ) and (15). More interesting is the case in which the four- 



fermion operator in Eq. < J13[ ) is generated at tree level, while the magnetic-moment transition in Eq. @ 
is generated only at one loop, as in models with imparity violation [14] or with leptoquarks p3[]. In this 
case, we expect 
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If Eq. ([16]) holds and if we take Mp — 1 TeV, then the ratios in Eqs. (|14|) and (p3|) become of order unity, 
so the different rare muon processes have comparable rates. 



Alternatively, if the fermion / in Eq. (13) is an electron (or a light quark), the effective operator 
can mediate \i — > 3e (or \i~-e~ conversion) at tree-level, and the coiTesponding process can dominate 



over the others [29]. For instance, we obtain 
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for the case / 

Figure [j] summarizes the behaviour of the ratio of branching ratios as a function of the relative 
strength of the effective operators in Eq. (||) and Eq. (|l3|), when the fermion / in Eq. (13) is an electron, 
as in the left part of Fig. [j], or a combination of first generation quarks, as in the right part of Fig. 111. It can 
easily be seen from the plots that when the magnetic moment operator dominates (A 2 <C A F ) the ratio 
of branching ratios saturates at several times 10 -3 , while it grows like (A 2 /A|,) 2 when the four-fermion 
operators are dominant (A 2 3> A F ). Interference effects are largest when A 2 ~ A 2 F , as expected. 




Fig. 1: Branching ratios normalised to — > ej) as a function of the ratio of the couplings of effective 
dimension-5 and dimension-6 operators (see text), for /i + — > e + e~e + (left) and /i~ N — > e~iV conversion in 
48 Ti. The solid (dashed) curves apply when the two operators interfere constructively (destructively). 



In conclusion, the various rare muon processes are all potentially very interesting. In the event of a 
positive experimental signal for muon-number violation, a comparison between searches in the different 
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channels and the use of the effective-operator approach discussed here will allow us to quickly identify 
the correct class of models. 

The study of rare muon decays also allows for the possibility of observing 'weak' time-reversal 
invariance violation (T violation)^]. While one cannot construct experimentally accessible T-violating 



observables in fi + — > e + 7 decay, in /i + — ► e + e~e + this is a possibility if the muon is polarized [30, 21]. 
In this case, one may test whether the angular distribution of the decay electrons depends on the T-odd 
Lorentz invariant oc ■ (p e ± x p e 2) (here p e j, i = 1, 2 are the momenta of the two positrons in the 
fi + rest frame and P„ is the muon polarisation). In particular, a T-odd asymmetry can be denned (see 



Refs. [30, 31] for details). In a number of theoretical models (see Refs. [31, 29, 32] and references 



therein) 'large' values (up to 24% [29]) for this T-odd asymmetry can be obtained. Unfortunately, in 
these cases the branching ratio for fi + — ► e + 7 is significantly larger than that for /i + — ► e + e~e + . 

It should also be noted that, again when the decaying muons are polarised, a number of P-odd 
observables in fi + — ► e + 7 and fi + — > e + e~e + can be defined. The measurement of these P-odd asym- 
metries can play a significant role when it comes to distinguishing different classes of models that yield 



significant flavour violation in the charged lepton sector (see Refs. [ 3jJ , [29J, [32|] and references therein). 



3.1.2 Rare muon processes in supersymmetry 

Supersymmetric extensions of the SM provide a very promising way of rendering the hierarchy of phys- 
ical mass scales more natural. Moreover, low-energy supersymmetry often leads to large sources of 
individual lepton number violation. It provides a framework for computing physics observables in a con- 
trolled manner as a function of a well-defined set of parameters. Those, in turn, can thus be constrained 
by the experimental limits. Whereas lepton-flavour number violation in rare \i decays may well have the 
same source as neutrino oscillations, the rates generically are no longer suppressed by powers of neu- 
trino masses: It is worth recalling that, for large values of the ratio of the two Higgs vacuum expectation 
values tan 0, the coefficient of the non-diagonal magnetic-moment operator in Eq. @ grows linearly 
with tan since, at the loop level, the chiral flip can be generated by the Higgs with the larger vacuum 
expectation value, and thus 

B(jjl -> erf) oc tan 2 /?. (18) 
We distinguish two different supersymmetric sources of lepton-flavour number violation. 
• Flavour non-diagonal soft terms 



It is well known [33] that a mismatch in flavour space between the lepton and slepton mass ma- 
trices generates tree-level transitions between different leptonic generations, both in charged and neutral 
currents. For instance, if we indicate the mixing angle between the first two generations of sleptons 
schematically by 9^, we obtain 

a 3 7T0f fi 

B( ^^c#- <19) 

where rh is a typical supersymmetric mass. In general, Og/i receives contributions from the flavour mis- 
matches of left- and right-handed sleptons and sneutrinos. Complete formulae for the rates of rare muon 
processes with the functional dependence on the different supersymmetric parameters can be found in 
Ref. [34]. As is apparent from Eq. (|l9|), the rates for rare muon processes can be large in supersymmetric 
models, since 0g/i is not necessarily suppressed by powers of neutrino masses. 

For generic values of the soft supersymmetry-breaking mass scale m ~ 1 TeV, the experimental 
limits on rare muon decays transform in very stringent upper limits on Ogp,. Therefore, it is common 

'it should be kept in mind that the study of CP violation is not practical due to the experimental difficulty of stopping 
negative muons in matter without absorbing them. 
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to invoke some universality condition or flavour symmetry that implies O^jx = at the scale at which 
the soft terms are generated. If this scale is sufficiently low, as in gauge-mediated supersymmetry, then 
loop corrections are small and rare muon processes are unobservable [|3~5|]. However, if the scale of 
supersymmetry breaking is large, as in supergravity models, then at the quantum level the soft terms 
are corrected by any high-energy flavour-violating interactions, and we generally expect more sizeable 



slepton mixing angles [ ]36| ] . 

For quantitative estimates we consider two examples for flavour-violating interactions at high 
energy. Both are commonly present in supergravity models. Additional sources of flavour violations 
should be expected and therefore our estimates should be viewed as conservative lower bounds. 

The first mechanism comes from renormalisable Dirac neutrino Yukawa couplings y v at energies 



larger than Mr [34, 37, 38, 39], the scale of the right-handed neutrino masses. Here we are assuming 



that the smallness of the neutrino masses results from the seesaw mechanism Q40Q . In the basis in which 
the charged-lepton Yukawa coupling matrix is diagonal, loop corrections to the slepton mass matrix are 
not diagonal in flavour space: 



87r 2 \v»)ki \y»)kj m Mr 



(20) 



Here mo and Aq are the universal soft supersymmetry-breaking scalar mass and trilinear term, respec- 
tively, and we have assumed that all three right-handed neutrinos have the same mass M r. The contri- 



bution in Eq. (20) induces flavour- violating mixing angles for both sneutrinos and charged sleptons, and 
therefore rare muon processes are generated by loop diagrams involving charginos and neutralinos. 



MSSM with right-handed neutrinos 
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Fig. 2: B(/j, — > e~f) as a function of the left-handed se- 
lectron mass for SUSY models with heavy right-handed 
neutrinos. The assumption here is that the right-handed 
neutrino mass matrix is proportional to the identity ma- 
trix with Mr = 10 13 GeV, and that all the neutrino mix- 
ing comes from the Yukawa couplings. LMA, SMA, and 
LOW refer to different solutions to the solar neutrino puz- 
zle (see Ref. @). 



The experimental information on neutrino oscillation parameters []41[], is not sufficient to recon- 
struct the complete structure of the relevant mass matrices. Therefore, it is necessary to rely on specific 
ansatze for them. Figure |2| depicts the value of B{[x — > ej) as a function of the left-handed selec- 
tron mass, assuming that the right-handed neutrino mass matrix is proportional to the identity and that 
Mr = 10 13 GeV. Results are shown for three different solutions to the solar neutrino puzzle and different 
values of the U e s element of the neutrino mixing matrix (see Ref. pi]]), assuming maximal *-> v T mix- 
ing as the solution to the atmospheric neutrino puzzle. Although the prediction for B(/j, — > ej) can vary 
significantly, depending on differences in the assumed texture of neutrino masses, supergravity models 
in which neutrino masses are obtained via the seesaw mechanism can induce fi — > at rates close to 
the current experimental sensitivity. This is especially true if the solution to the solar neutrino puzzle is 
indeed the large mixing angle solution or if the U e 3 element is nonzero. It should be noted, however, 
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that because the magnitudes of neutrino Dirac Yukawa couplings are related to the right-handed neu- 
trino mass scale Mr in seesaw models, B(/j, — > cy) is very sensitive to the right-handed neutrino scale, 
roughly 

B(n -> erf) oc M R [ln(M Pl /M R )} 2 . (21) 

As can be seen in Fig. ^, the prediction for B(/j — ► ej) also depends quite sensitively on the assumed 
value of U e 3. This is a reflection of the general sensitivity to assumptions on the textures of the charged- 
lepton and neutrino mass matrices. 
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Fig. 3: Scatter plot of values of — ► e7) and 
B(t — ► ^7) in a sampling of models of lepton masses 
and mixing with textures generated by £7(1) flavour sym- 
metries. There is considerable spread in the predictions, 
but it seems that fj, — > e7 decay may offer better prospects 
for discovering new physics than does r — > ^7, within 
this class of models p^]. Note that current experiments 
already constrain — * 67) < 1.2 x 10" 11 and 



Figure |3| displays values of B(fi — > cy) and B(t — > ^7) found [39] in a sampling of specific 
models whose masses and mixing angles have different textures generated by £7(1) flavour symmetries. 
We see that the predictions span a considerable range, but, in this class of models, values of B(fi — > erf) 
close to the present experimental limit tend to be associated with values of B(t — > fij) that are con- 
siderably below the present experimental upper limit B(t — > ^7) ~ 10~ 6 , suggesting that B(n — > e^f) 
may be a more promising place to search for new physics. On the other hand, improving the sensitivity 
to B(t — > ^7) provides very useful information when it comes to disentangling different neutrino mass 



models [38] 



The decay — ► e + e e + is usually dominated by an intermediate photon, and therefore Eq. (11) 



holds rather generally. This dominance is less certain in fi~-e~ conversion, so there may be some 
deviations from Eq. (fl^). One particular example [[^] is shown in Fig. |], where the penguin-diagram 
contribution (which includes the magnetic-moment part) is plotted separately from the sum of penguin 
and box diagrams. As in Fig. [I[ we see that Eq. ( fl2| ) holds generically to within some numerical factor 
of order unity. 




Fig. 4: The ratio of the rates for /i-e conversion and 
H — > e-f, as a function of mo for \i > and /i < 
(tan/3 = 3) [^]. The penguin contribution, which in- 
cludes the magnetic -moment piece, is plotted separately 
from the sum of the penguin and box diagrams. The re- 
sults are in the range expected from Eq. (n2h. 



The possible discrepancy between the Standard Model and the recent measurement of the muon 
anomalous magnetic moment 5a^ = (g^ — 2) [ |42| ] puts a new perspective on the searches for fj, — > e^f 
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and related processes. Quite generally, the discrepancy, if confirmed, would be evidence for new physics 



at the TeV scale that might be detectable in other processes involving muons [}43|]. Within any specific 
model that incorporates lep ton-flavour violation as inferred from neutrino-oscillation experiments, the 
magnitude of 5a^ may be used to normalise the prediction for B(/i — > erf). This correlation is very close 
if 8a n is due to supersymmetry [|4|, 45, 16], since fi — > ej decay is dominated by the magnetic-dipole- 
moment type operator in Eq. @, and the supersymmetric contribution <5a^ USY to the muon anomalous 
magnetic moment is induced by similar diagrams mediated by sleptons, neutralinos and charginos. If 
the left-handed sleptons have lepton-flavour-violating soft masses, the correlation is particularly strong, 
since the dominant contributions to both H^erf and daf^ are from the left-handed slepton-chargino 
loop diagrams: 



SUSY 



5a 



5a| + a\ 



48vr 



??? 



tan/?, 



SUSY 



7T 5 

— a(a 2 + ^ay) 2 (G|m| US Y) _1 tan 2 f3 

, fSa^V A™| )i2\ 2 
3 x 10~ 5 ' ^ 1 ' tL 



m 



SUSY 



I 10 -° ) 



m 



(22) 
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SUSY 



Here we have assumed, for illustration, that all the supersymmetric mass parameters are equal to a com- 
mon mass, mgusY- Assuming that the 2.6-sigma deviation recently observed in the Brookhaven E821 
experiment | |42| ] is due to low energy SUSY, (i.e.5a^ SY ~ 10~ 9 ), we can use the present experimental 
upper limit B(fi — ► &y) < 1.2 x 10 -11 to obtain a stringent constraint on the lepton-flavour-violating 
mass insertion: 
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(25) 



As can be seen from Eq. (20), for example, large neutrino Yukawa couplings could easily induce an 
observable effect in the near future. The correlation between B(fi — > ef) and <5a^ USY in the case of 
a specific choice of neutrino mass texture p5j], is illustrated in Fig. ||. We see that, at least within the 
context of this model, [i — > e-f decay may well be observable in the next round of experiments, and the 
same is true for anomalous /x — ► e conversion on nuclei. The establishment of the anomaly in — 2, 
together with improvements of the present experimental limits on rare muon decay processes, will probe 
severely supersymmetric models. 
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Fig. 5: The contours B(n -> e<y) = 10 -11 , 10" 12 , KT 13 and 
10 -14 , for a specific neutrino mass texture described in Ref. [|45|l, 
in the gaugino (myj) versus scalar mass (mo) plane, assuming 
tan/3 = 10, A = 0, and fj, > 0. The regions allowed by the 
E821 measurement of a M at the 2-a level [ p^ ] are shaded (pink) and 
bounded by solid black lines, with dashed lines indicating the l-er 
ranges. The dark (brick-red) shaded regions are excluded because the 
LSP is the charged stau and the light (turquoise) shaded regions are 
those with 0.1 < fl x h 2 < 0.3 that are preferred by cosmology. Also 
shown are the contours nth = 113 GeV (assuming m t = 175 GeV) 
andm x + = 104 GeV. 



The second mechanism we investigate is present when the minimal supersymmetric extension 
of the SM is embedded into a grand unified theory (GUT). In this case, leptons and quarks belong to 
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the same group representation, and the large top Yukawa coupling y t will split the mass of the third- 
generation slepton from the other two [ 47| , 48]. In SU(5) models [|^, 50], the up-type quarks are in the 
same representation as the right-handed charged leptons, and therefore we find 



777; 



33 



3(3777.2 +.4.2) Mpi 

-Vt m ■ 



Mqut 



(26) 



Once the charged-lepton mass matrix is diagonalised, Eq. (|26j) induces slepton mixing. The left part of 
Fig. H depicts B(fi — > erf) as a function of the right-handed selectron mass for different values of tan (3 in 
the case of an SU(5) GUT. In this case, since the mixing is present only in the right-handed lepton sector, 
the positron in the anti-muon decay is left-handed, /j + — > e^j. This decay could be distinguished from 
/j + — > e^7 by studying the angular distribution of the positron in polarised /i + — > e + j decays [^|, [52| ]. 
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Fig. 6: B(fj, — > ef) as a function of the right-handed selectron mass in supersymmetric models with SU(5) (left) 
M\ and SO(10) (right) @ grand unification. 



In the case of SU(5) models with massless neutrinos, the rates for fi — ► e^f are typically much 
smaller than the present experimental bound. This is because chargino loop diagrams are absent and 
neutralino diagrams suffer from a partial cancellation. Note that the position of the dip in the curves of 
Fig. ^ (corresponding to the exact cancellation) depends on the value of the gaugino mass M\. However, 
if higher-dimensional operators are added to models with SU(5) symmetry, B(fj, — > cy) may be enhanced 
(see, for example, Ref. [|50|]). Such operators may be required in order to obtain realistic mass relations 
for first and second generation quarks and leptons. Also note that the values of B{fjL — ► &y) depends very 
strongly on the precise value of the top Yukawa coupling at the GUT scale, which is not well determined 
by the measurement of the top quark mass, because of its quasi-fixed-point structure. 

The situation changes in SO(10) models^ [ |48| , 54], where all SM fermions belong to the same 
irreducible representation, and therefore the top Yukawa effect splits not only £r, but also £l- Further- 
more, since the loop diagram involves the exchange of a third-generation slepton, the necessary chiral 
flip can occur along the f line, and B(pL — ► erj) is enhanced by a factor m^/m 2 ^. Figure fright) depicts 
B(fj, — > ej) as a function of the right-handed selectron mass for different values of tan @ in the case of 
anSO(10) GUT. 



Other unification scenarios have also been explored. See, for example, Ref. 
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• R-parity violation 

In supersymmetric models where R parity is not imposed as a discrete symmetry, one encounters 
renormalisable interactions that violate total lepton number and individual lepton numbers. At the scale 
of the muon mass, these interactions correspond to effective four-fermion operators of the kind discussed 
in Section 3.1.1. Therefore, depending on the field content of the operator, we can have situations in 
which /i + — ► e + e~e + or fi~-e~ conversion become the dominant rare muon process. For instance, if a 
model predicts the simultaneous presence of the two i?-parity- violating interactions X\3iL^/ L^' E^ 1 ' 
and XmiL^ , then fi + — ► e + e~e + occurs at the tree level with a branching ratio 

B(» - 3e) = = ( 2"MV 10- . (27) 

Present constraints from rare muon processes provide very stringent upper limits on numerous products 
of i?-parity-violating couplings, some of which are listed in Table [j]. At a neutrino factory complex, 
these bounds could be improved by two orders of magnitude if fj, — > ej and fi — > eee searches are not 
successful and by three orders of magnitude if fi-e conversion in nuclei is not observed. 

Table 1: Upper limits on products of .R-parity- violating couplings from the current limits for muon-number- 
violating processes, taking slepton masses equal to 100 GeV and squark masses equal to 300 GeV. The values 
in parentheses indicate the sensitivity which could be achieved at a neutrino factory complex, i.e., resulting from 
improvements in experimental sensitivities to rare fi decays by 4-6 orders of magnitude. The coupling constants A 
and A' refer to the interactions D XijkL l L J E k + X'^ISQi D k , where i,j,k are generation indices. Tree-level 
constraints are indicated by [tree]. 
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3.1.3 Rare muon processes in models with extra dimensions 

Recently, theories with extra spatial dimensions have attracted a lot of attention as a possible solution 
to the problem of the large hierarchy between the Planck and Fermi mass scales. The hypothesis is that 
the universe possesses 1 + n (n > 3) space-time dimensions, while the Standard Model particles are 
constrained to live on a 1+3-dimensional subspace. Gravity, which is described by the geometry and 
therefore propagates in the full space, appears very weak to us either because its strength is diluted in a 
large compactified extra-dimensional space [55] or because it is localised away from us in spaces with 
non-factorizable geometries [pq]. Since these theories assume that gravity becomes strongly coupled 
at an energy scale Mq comparable to the electroweak scale, the possibility of explaining the smallness 
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of the neutrino masses via the classical see-saw mechanism is lost. Nevertheless, the small neutrino 
masses could now have an explanation based on geometrical arguments, similar to those that led to a 
justification of the small ratio Mw/Mpi. For this to happen, one needs to assume the existence of right- 
handed neutrinos which, like the graviton, also propagate in the extra dimensions. If this is the case, their 
Yukawa interactions with the SM left-handed neutrinos are effectively suppressed by large geometrical 



factors []57[], and therefore the SM neutrinos obtain, after electroweak symmetry breaking, a very small 
Dirac mass. When neutrino family mixing is included, one finds that the Kaluza-Klein modes of the 
right-handed neutrino mediate, at the loop level, flavour transitions in the charged sector [p 



,59] 



Particularly interesting are 'minimal models' in which all flavour transitions are described by only 
two free parameters, besides the observable neutrino masses and mixing angles. These two parameters 
are the 'fundamental' cut-off scale A (expected to be of the order of the weak scale, if these models are 
motivated by the hierarchy problem) and a dimensionless coefficient e, which is currently constrained 
to be approximately less than 10~ 2 and whose expected value depends strongly on the neutrino Yukawa 
couplings and on details of the extra dimensional model, such as the number of extra dimensions 6 in 
which the right-handed neutrinos propagate (see Ref. [ ]59| ] for details). 

Under these conditions, the branching ratio for \i 



e7 is 



J e2U*2-7 rp 2 - + U e3 U* 3 

p Am 



atm 



(28) 



Here the unitary matrix U describes the neutrino mixing angles and Am 2 un atm are the neutrino mass- 
squared differences relevant to solar and atmospheric oscillation, for a hierarchical neutrino mass spec- 
trum. 

Unlike ti — ► £jj decays, the rates for ii — ► eee and \i — > e conversion in nuclei are quite dependent 
on the unknown ultraviolet details of the models. Nevertheless, these rates can be predicted as a function 
of e and A (see Refs. [E 



59] for complete expressions). It should be noted that the rates for li 



eee and fj, — > e conversion in nuclei can be significantly enhanced with respect to the rate for fx — > 
e7 in some regions of the e, A parameter space, similarly to the case of SUSY models with imparity 
violation. Figure [7] shows the rates for the different flavour-violating lepton processes, as functions 
of |f7 e3 |, assuming Am s 2 m /Am^ tm = 1CT 2 , \U^/U r3 \ 2 = 1, \U e2 /U e i\ 2 = 2/3 (as suggested by 
atmospheric and LMA solar oscillations), and no CP violation in the neutrino mixing matrix, for different 



values of e and A. We also show the conversion probabilities for ia. 



; e and 



v T at very short 



baseline experiments. This allows a comparison of the discovery reach for different facets of a neutrino 
factory, namely, experiments with stopped muons and experiments with an intense neutrino beam. 

We emphasize that the one-loop effects giving rise to rare muon and tau processes that violate indi- 
vidual lepton flavour number are cut-off dependent and can only be qualitatively estimated. Nonetheless, 
in minimal models, their flavour structure is directly related to the physical neutrino mass matrix, and 
one is capable of predicting and relating the rates of rare charged and neutral lepton processes in terms 
of observable neutrino oscillation parameters. This is in sharp contrast to other cases of physics beyond 
the SM. In SUSY models, for instance, the rates for rare muon processes are perturbatively calculable, 
but their relations with neutrino oscillations parameters are indirect and strongly model-dependent, as 
described in the previous subsection. 
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Fig. 7: Branching ratios of /i — > e7, fi — > 3e and r — > ^7, normalized rate of /x-e conversion in 27 A1, and 
transition probabilities for 1^ — > v T and f M — ► v e at very short baseline as functions of [ C/e3 1 . We have chosen 
the parameters 5 = 5, e = 0.003 (top) or 0.0003 (bottom), and A = 1 TeV (left) or 10 TeV (right). We assumed 
Am s 2 un /Amj tm = 10 -2 , IU^/Uts] 2 = 1, \U e 2/U e \\ 2 = 2/3 (i.e. maximal mixing in the atmospheric sector and 
the LMA solution to the solar neutrino puzzle), no CP-violation in the neutrino mixing matrix, and hierarchical 
neutrino masses (m 2 <C m| <C m 2 ). See Ref. J59| ] for details. 
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Searches for Lepton Number Violation 
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Fig. 8: Historical development of the 90% C.L. upper limits (UL) on branching ratios respectively conver- 
sion probabilities of muon-number violating processes which involve muons and kaons. Sensitivities expected 
for planned searches are indicated in the year 2008 (see also Ref. p^]). The projections for a neutrino factory 
(NUFACT) are also shown. 



3.2 Experimental prospects 

The experimental sensitivities achieved during the past decades in tests of muon number conservation 
are illustrated in Fig. |8[ Generally the tests were limited up to now by the intensities of the available fi 
and K beams, but in some cases detector limitations have played a role as well. 



All recent results with /i + beams were obtained with 'surface' muon beams [p0|], i.e. beams of 
muons originating in the decay of 7r + 's that stopped in the pion production target. These beams offer 
the highest muon stop densities that can be obtained at present, allowing for the low-mass experimental 
targets that are required for the ultimate resolution in positron momentum and emission angle or the 
efficient production of muonium in vacuum. 

In this subsection we study the question of how far experimental searches could benefit from 
muon beam intensities which are 2-3 orders of magnitude higher than presently available. We do so by 
analysing the present state of the art, assuming modest improvements in detector technology. 

3.2.1 fj, — * e7 

During the past 25 years the experimental sensitivity to this decay mode was raised by two orders of 
magnitude (see Table §). The most sensitive search was performed recently by the MEGA Collaboration, 



Table 2: 90% C.L. upper limit on the branching ratio for ji — > obtained by previous experiments at meson 
factories 



Experiment or Lab. 


Year 


Limit 




Ref. 


SIN (presently PSI) 


1977 


1 x 10" 


y 




TRIUMF 


1977 


3.6 x 10" 


9 




LANL 


1979 


1.7 x 10" 


-10 




Crystal Box 


1986 


4.9 x 10' 


-11 




MEGA 


1999 


1.2 x 10' 


-11 


§1 
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establishing an upper (90% C.L.) limit on B^_ >C7 of 1.2 x 10 ||65|]. An experiment [|66|], which aims 
at a single-event sensitivity of ~ 10~ 14 was approved by the PSI scientific committee in July 1999. 
Measurements should start in 2003 and data taking should go on for two or three years. This experiment 
will use a surface muon beam that reaches an intensity around 5 x 10 8 fi + / s. As a next step it seems 
reasonable to consider experiments aiming at a sensitivity of 10~ 15 or better. As we shall see below, 
it is not at all obvious how to reach such levels of sensitivity without running into the background of 
accidental e^j coincidences. 
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Experimental sensitivity 

In order to understand how the sensitivity^] to the \i — > e 7 decay is related to the experimental 
resolutions, let us consider an experiment in which muons are stopped in a thin (typically 10 mg/cm 2 ) 
target; photons and positrons are detected by a pair spectrometer or an electromagnetic calorimeter and 
a magnetic spectrometer, respectively The expected number N s of observed n — ► e 7 decays can be 
written as: 

n 

N s = R^T— e e e y e cut B^_» e7 , (29) 

where is the muon stop rate, T is the total measuring time, Q is the detector solid angle (we assume 
identical values for the photon and the positron detectors), e e and e 7 are the positron and photon detection 
efficiencies, e cu t is the efficiency of the selection cuts. Selection cuts can be applied on the reconstructed 
positron energy (E e ), photon energy ( Ej ), opening angle (9 ei ) and relative timing (t ej ). Considering 
selection windows (indicated with a A in front of the relevant quantity) covering 90 % of the signal (full 
window = 1.4 FWHM for Gaussian distributions) we have e cu t = 0.9 ~ 0.66. 

The level of prompt background from fi + — > e + vv^ depends on the selection windows on the 
energies and opening angle of the e^j pair. The background level predicted with no angular cut is shown 
in Fig. ^. Whereas these estimates do not take into account details of the detector response functions, 
they are accurate enough for our purpose. It appears that detector technology is sufficiently advanced to 
keep prompt background at a negligible level 




Fig. 9: Prompt background from /i + — > e + vv^ as a function of the selection windows on the positron (x = 
2E e + /m^) and photon (y = 2E~ f /m tl ) energies. We have defined 8x = 1 — x and Sy = 1 — y (results from 
Ref. H). 

In the proposed PSI experiment (with ~ 10 8 /u/s) the background is dominated by accidental 
coincidences of a positron from the normal muon decay and a photon which may originate in the decay 
/i + — > e + vV^ or may be produced by an e + through external bremsstrahlung or annihilation in flight. In 
a DC beam the number of accidental coincidences is given by: 

N b = R^f e e e f^(— ) 2 7t-^l2 At T , (30) 

where f e (/ 7 ) is the e + (7) yield per stopped muon within the selection window and At is the cut applied 
on the e + — 7 time difference. For a non DC beam N(, must multiplied by the inverse of the duty cycle. 

3 In this document, sensitivity is defined as the 90% CL upper limit in the absence of a signal. 
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An evaluation of f e and / 7 is given in Refs. Q67| , |68[ |. Figure [10| shows / 7 as a function of the 
energy threshold calculated for a 50mg/cm 2 target. 




Fig. 10: Photon yield vs. energy threshold. Dotted: annihilation in flight in a 50 mg/cm 2 target. Dashed: 
radiative decay. Solid: sum of the two. 



By taking into account the Michel spectrum for the positron and only the radiative muon decay for 
the photon, one finds near the endpoint (E e and£^ ~ m^/2): 
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From Eqs. (|29|) and ( J3QD the signal-to-background ratio is 

^cut-B^t— > e 7 
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(32) 



The detector resolutions quoted in the PSI proposal are shown in Table |3| together with those of 
previous experiments. For this experiment a solid angle Si/An of about 10% and a total measuring time 



Table 3: Detector resolutions, muon stop rates and accelerator duty cycles of the various /.t — > e 7 searches. 



Exp. /Lab. 


AE e /E e 


AEy j Ery 






Instantaneous stop rate 


Duty cycle 




% FWHM 


% FWHM 


(ns) 


(mrad) 


(1/s) 


(%) 


SIN (now PSI) 


8.7 


9.3 


1.4 




(4 ... 6) x 10 5 


100 


TRIUMF 


10 


8.7 


6.7 




2 x 10 5 


100 


LANL 


8.8 


8 


1.9 


37 


2.4 x 10 6 


6.4 


Crystal Box 


8 


8 


1.3 


87 


4 x 10 5 


(6. ..9) 


MEGA 


1.2 


4.5 


1.6 


17 


2.5 x 10 s 


(6. ..7) 


PSI 


0.7 


1.4 


0.15 


12 


10 8 


100 



T of almost a year are being considered. It turns out that one background event is expected and one 
fi — ► 6 7 event for B^ ei = 10 -14 . 

Since the accidental background rises quadratically with the muon stop rate, it will be even more 
problematic in future experiments using a higher beam intensity. An experiment with = 10 10 /u/s 
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and all the other quantities of Eq. (29) unchanged would yield one p — > event for B^ ej = 10 16 . 



However, the accidental background would increase to 10 4 events. It is obvious that better detector 
resolutions so as to reduce the product 

^(^AO^At (33) 

e "y 

and/or improved experimental concepts are required. 
Ideas for background reduction 

We report here on some possibilities discussed in the study group. 
Detector improvements 

The safest improvement can be obtained in the positron momentum resolution where a substantial 
factor (~ 10) can be gained by utilising a beta-spectrometer [ |67| ] which may reach resolutions of the 
order of 0.1%. As an additional advantage such a spectrometer would result in reduced rates in the 
tracking detectors as compared to the broad-band systems used so far. 

The resolution in e + — 7 opening angle is determined by positron multiple scattering in the stop- 
ping target, so another possible improvement could be the reduction of the target thickness. The use 
of very intense muon beams would give the possibility to obtain a 10 10 /x/s stop rate with much thinner 
targets. A reduction by one order of magnitude would reduce A9 ey by a factor of 3, corresponding to 



an order of magnitude less accidental background, see Eq. (|33|). Photon background from positron anni- 
hilation would be reduced by the same factor. Since the decay fi + — > e + vl7j would then dominate the 
photon background, one could envisage reducing even this background by vetoing on the accompanying 
positron. 

Target subdivision 

Another way to diminish the accidental background could be the use of a row of very thin targets. 
The requirement for positrons and photons to originate in the same target would reduce the background 
in proportion to the number of targets. One might think of it as a series of identical p — > e 7 set-ups 
sharing the same beam. However, this configuration poses several complications: 

• the photon detector must have sufficient directional selectivity to distinguish between two adjacent 
targets; 

• it is not obvious how to combine the scheme with the idea of a high-resolution positron spectrom- 
eter. 



Muon polarization 

The use of the muon polarization for reducing the accidental background has been suggested in 
the past [|| In muon decay, positrons and photons at energies close to the endpoint follow a 
(1 + P^ cos 9) angular distribution (where is the degree of muon polarization). Configurations with 
back-to-back e + -7 pairs are therefore suppressed with respect to the unpolarized case. If one requires 
a reasonable solid angle subtended by the detectors (Q/Air ~ 0.1 ) the rejection factor turns out to be 
approximately 5. It should be noted, however, that in some models \i — > e 7 would be suppressed as 
well. 
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3.2.2 fi+ -> e+e+e" 

From an experimental point of view the decay \i — > 3e offers some important advantages compared to 
the more familiar ll — ► discussed in the previous section. Since the final state contains only charged 
particles there is no need for an electromagnetic calorimeter with its limited performance in terms of 
energy and directional resolution, rate capability, and event definition in general. Apart from the obvious 
constraints on relative timing and total energy and momentum, as can be applied in ll — > e^f searches as 
well, there are powerful additional constraints on vertex quality and location to suppress the accidental 
background. Of major concern, however, are the high rates in the tracking system which has to stand the 
load of the full muon decay spectrum. 

The present experimental limit, B(fi — ► 3e) < 1 x 1CP 12 [|6^], was published in 1988. Since 
no new proposals exist for this decay mode we shall analyse the prospects of an improved experiment 
with this SINDRUM experiment as a point of reference. A detailed description of the experiment may 
be found in Ref. 0. 

Data were taken during a period of six months in the years 1984 and 1986 when a 25 MeV/c 
subsurface beam was brought to rest in a hollow double-cone target at a rate of 6 x 10 6 ll + s _1 . The 
target geometry is described in Table Q The spectrometer acceptance for ll — > 3e was 24% of 4ir sr (for 



Table 4: SINDRUM I target for ll -> 3e 



Shape 


hollow double-cone 


Material 


foam 


Length 


220mm 


Diameter 


58mm 


Mass/cm 2 


llmg/cm 2 


Total mass 


2.4 g 



a constant transition-matrix element) so the only place for a significant improvement in sensitivity would 
be the beam intensity. SINDRUM I is a solenoidal spectrometer with a relatively low magnetic field of 
0.33 T corresponding to a transverse-momentum threshold around 18 MeV/c for particles crossing the 
tracking system. This system consists of five cylindrical MWPCs concentric with the beam axis. Three- 
dimensional space points are found by measuring the charges induced on cathode strips oriented ±45° 
relative to the sense wires. Gating times were typically 50 ns. 



Figure |11] shows the time distribution of the recorded e + e + e~ triples. Apart from a prompt contri- 
bution of correlated triples one notices a dominant contribution from accidental coincidences involving 
low-invariant-mass e + e~ pairs. These are explained by Bhabha scattering of positrons from normal 
muon decay ll — ► evV. The accidental background thus scales with the target mass, but it is not obvious 
how to reduce this mass significantly from the 11 mg/cm 2 shown in Table Q 

Figure 12 shows the vertex distribution of prompt events. One should keep in mind that most of 
the uncorrelated triples contain e + e~ pairs coming from the target and their vertex distribution will thus 
follow the target contour as well. This 1-fold accidental background is suppressed by the ratio of the 
vertex resolution (couple of mm 2 ?) and the target area. There is no reason, other than the cost of the 
detection system, not to choose a much larger target than that given in Table |j. Such an increase might 
also help to reduce the load on the tracking detectors. 

For completeness, Fig. 13 shows the distribution of total momentum versus total energy for three 
classes of events, (i) uncorrelated e + e + e~ triples, (ii) correlated e + e + e~ triples, and (iii) simulated 
ll — > 3e decays. The distinction between uncorrelated and correlated triples has been made on the basis 
of relative timing and vertex as discussed above. 

What would a ll — > 3e set-up look like that would aim at a single-event sensitivity around 10~ 16 , 
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high 



Fig. 11: Relative timing of e 



events. The two positrons are labelled according to the invariant mass when 



combined with the electron. One notices a contribution of correlated triples in the centre of the distribution. These 
events are mainly // — > Zevv decays. The concentration of events along the diagonal is due to low-invariant-mass 
e + e~ pairs in accidental coincidence with a positron originating in the decay of a second muon. The e + e~ pairs 
are predominantly due to Bhabha scattering in the target. 

i.e., would make use of a beam rate around 10 10 /x + /s? The SINDRUM I measurement was background- 
free at the level of 10~ 12 with a beam of 0.6 x 10 7 fi + /s. Taking into account that background would 
have set in at 10~ 13 , the increased stop rate would raise the background level to « 10~ 10 , so six orders of 
magnitude in background reduction would have to be achieved. Increasing the target size and improving 
the tracking resolution should bring two orders of magnitude from the vertex requirement alone. Since 
the dominant sources of background are accidental coincidences between two decay positrons (one of 
which undergoes Bhabha scattering) the background rate scales with the momentum resolution squared. 
Assuming an improvement by one order of magnitude, i.e., from the « 10% FWHM obtained by SIN- 
DRUM I to ~ 1% for a new search, one would gain two orders of magnitude from the constraint on total 
energy alone. The remaining factor 100 would result from the test on the collinearity of the e + and the 
e + e~ pair. 

As mentioned in Ref . J70| ] a dramatic suppression of background could be achieved by requiring 
a minimal opening angle (typically 30°) for both e + e~ combinations. Depending on the mechanism for 
fi — > 3e, such a cut might, however, lead to a strong loss in \i — > 3e sensitivity as well. 

Whereas background levels may be under control, the question remains whether detector concepts 
can be developed that work at the high beam rates proposed. A large modularity will be required to solve 
problems of pattern recognition. Also the trigger for data readout may turn out to be a great challenge. 



3.2.3 fie conversion 

Neutrinoless \x~ — > e~ conversion in muonic atoms, fi~(A, Z) — > e~(A, Z) with A mass number and Z 
atomic number, offers some of the best tests of lepton flavour conservation (LFC). For conversions leav- 
ing the nucleus in its ground state the nucleons act coherently, which boosts the conversion probability 
relative to the rate of the dominant process of nuclear muon capture. 

Assuming LFC violation one may wonder how the conversion probability B fie varies as a function 
of Z and A, and with what probability the nucleus stays in its ground state. Earlier calculations J7l| , f72| , 
73] predicted a steady rise of the branching ratio until Z 30, from where on it was expected to drop 
again. For this reason most experiments were performed on medium-heavy nuclei. More recently it has 
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Fig. 12: Spatial distribution of the vertex fitted to prompt e + e + e triples. One clearly notices the double-cone 
target. 




Fig. 13: Total momentum versus total energy for three event classes discussed in the text. The line shows the 
kinematic limit (within resolution) defined by T\pc\ + |Spc| < m M c 2 for any muon decay. The enhancement in 
the distribution of correlated triples below this limit is due to the decay fj, — ► 3ei/u. 



been estimated that may keep increasing with Z [74, 75]. The same calculations predict the coherent 
fraction to be larger than 80% for all nuclear systems. 

The dependence of B^ e on the normalized neutron excess (N — Z)/ (N + Z), with N = A — Z, 
depends on the nature of the LFC-violating propagator. Although no mechanism considered so far gives 
rise to a cancellation of the neutron and proton contributions for any given nucleus, a general model- 
independent analysis requires at least two measurements with significantly different values for (N — 
Z)/(N + Z)^,^]. 

When negative muons stop in matter, they quickly get captured and form muonic atoms, which 
mostly reach their ground state before decaying. The main decay channels are muon decay in orbit 
fi~(A,Z) — > e'VeV^A, Z) and nuclear muon capture fi~(A,Z) — ► v^A, Z — 1). Experimentally 
coherent n — > e conversion offers a number of advantages. Since the observation of the process requires 
the detection of only one particle in the final state, there are no problems with accidental coincidences, 
which constitute one of the major limitations in searches for fi — > and /i — ► 3e. The electron is 
emitted with a momentum p e m m^c, which coincides with the endpoint of muon decay in orbit (MIO), 
which is the only intrinsic background. Since the momentum distribution of muon decay in orbit falls 
steeply above the experimental set-up may have a large geometrical acceptance and the detectors 
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can still be protected against the vast majority of decay and capture events. Energy distributions for MIO 



electrons have been calculated for a number of muonic atoms | |77| ] and energy resolutions in the order of 
1% are sufficient to keep this background below 10~ 16 . 

There are several other potential sources of electron background in the energy region around 100 
MeV, involving either beam particles or cosmic rays. Beam-related background may originate from 
muons, pions or electrons in the beam. Apart from MIO, muons may produce background by muon decay 
in flight or radiative muon capture (RMC). Pions may produce background by radiative pion capture 
(RPC). Capture gammas from RMC and RPC produce electrons mostly through e + e~ pair production 
inside the target. 

Beam-related background can be suppressed by various methods: 

• Beam pulsing 

Since muonic atoms have lifetimes of order 100 ns, a pulsed beam with buckets short compared 
to this lifetime would allow one to remove prompt background by measuring in a delayed time 
window. As will be discussed below there are stringent requirements on the beam suppression 
during the measuring interval. This is the concept of the new MECO proposal. 

• Beam veto 

In the past, prompt background has been removed with the help of a beam counter. For the beam 
intensities proposed here this method can not be applied. 

• Beam purity 

A low-momentum (< 70 MeV/c) beam with extremely low pion contamination (< 100 7r~'s 
per day) would keep prompt background at a negligible level. A major advantage of the method 
is that heavy targets such as gold with lifetimes around 70 ns can be studied. This scheme was 
applied by SINDRUM II and will be discussed in the following. 



At present, fie conversion is being searched for by the SINDRUM II Collaboration at PSI y78Q 
and a new experiment (MECO [^]) is planned at BNL. The sensitivity levels are rs 5 x 10 -13 and 
w 5 x 10 -17 , respectively. 



Figure 14 shows the SINDRUM II spectrometer. Beam particles are extracted from a carbon 



production target at backward angles and momentum-analysed in a conventional quadrupole channel. 



The beam is focused into a 9 m long transport solenoid (see Fig. |14|) called PMC (pion muon converter). 

A pion reaching the gold target has a chance of order 10~ 5 to produce an electron in the energy 
region of interest, so the pion stop rate must be below one every ten minutes. At the PMC entrance 
the beam contains similar amounts of muons and pions. Since the pion range in matter is about half as 
large as the corresponding muon range, the pion contamination can be reduced strongly with the help 
of a moderator at the PMC entrance. Only one out of 10 6 pions may cross this moderator. Typically 
99.9% of them would decay before reaching the target. The requirement puts strong constraints on the 
high-momentum tail transported by the beam line which could be met after a careful optimisation of the 
beam settings. 

During an effective measuring period of 75 days about 4 x 10 13 muons stopped in the gold target. 



Figure [15] shows as a preliminary result various momentum distributions. The main spectrum, taken 
at 53 MeV/c, shows the steeply falling distribution expected from muon decay in orbit. Two events 
were found at higher momenta, but just outside the region of interest. The agreement between measured 
and simulated positron distributions from /i + decay gives confidence in the momentum calibration. At 
present there are no hints about the nature of the two high-momentum events: they might be induced by 
cosmic rays or RPC, for example. Both processes result in flat momentum distributions such as shown 
by the data taken at 63 MeV/c. 

As a preliminary result the single-event sensitivity was estimated to be » 2 x 10~ 13 , i.e., an 
improvement by two orders on magnitude of the previous best result on a heavy target 
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SINDRUM II 



Fig. 14: SINDRUM II 

The MECO experiment plans to combat beam-related background with the help of a pulsed 8 
GeV/c proton beam. Figure [R] shows the proposed layout. Pions are produced by 8 GeV/c protons 
crossing a 16 cm long tungsten target, and muons from their decays are collected efficiently with the help 
of a graded magnetic field. Negatively charged particles with 60-120 MeV/c momenta are transported by 
a curved solenoid to the experimental target. In the spectrometer magnet a graded field is applied as well. 
Figure [l7| shows the proposed time structure of the proton beam. A major challenge is the requirement for 
proton suppression in between the proton bursts. In order to keep the pion stop rate in the 'silent' interval 
below the aforementioned 100 per day, a beam extinction factor better than 10 _8 -10 -9 is required. 

Table |5| compares the \x~ intensity and single-event sensitivity obtained by SINDRUM II with 
those expected for MECO and the v factory. 

Table 5: Rates in fxe conversion searches 





SINDRUM II 


MECO 


NUFACT 




PSI 


BNL 


CERN 


Status 


Stopped 


Proposed 


Discussed 


Proton momentum [GeV/c] 


1.2 


8 


2 


Proton rate [s _1 ] 


1 x 10 16 


4 x 10 13 


1 x 10 16 


fi~ rate [s _1 ] 


3 x 10 7 


1.7 x 10 11 


10 14 


yT stops [s" 1 ] 


(1...2) x 10 7 


1.0 x 10 11 


2.0 x 10 13 


Single event sensitivity 


2 x 10~ 13 


2.0 x 10~ 17 


1Q -19 



Two scenarios for a \i beam are under study (see also Section 6.2 below): 
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Fig. 15: Recent results by SINDRUM EL Momen- 
tum distributions for three different beam momenta 
and polarities: (i) 53 MeV/c negative, optimised for 
fjT stops, (ii) 63 MeV/c negative, optimised for tt~ 
stops, and (iii) 48 Me V/c positive, for /i + stops. The 
63 Me V/c data were scaled to the different measur- 
ing times. The fi + data were taken at reduced spec- 
trometer field. 
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Fig. 16: MECO 



• A bunched proton beam extracted from the accumulator could be used in a fashion similar to the 
MECO beam line. As explained above, this option requires extinction factors of ten orders of 
magnitude or so to keep prompt background under control. 

• As an alternative one might use a cooled muon beam. Such a beam would probably be free of pions 
and contaminated only by decay electrons. As a second by-product of the cooling the extinction 
factor should be much higher at this stage. 



3.2.4 Muonium-antimuonium conversion 

Muonium is the atomic bound state of a positive muon and an electron. For leptons a spontaneous 
conversion of muonium (/i+e _ ) into antimuonium (/i~e + ) would be completely analogous to the well 
known K° — K° oscillations in the quark sector. A search already suggested in 1957 by Pontecorvo [81] 
three years before the atom was discovered by Hughes [|8"2j]. The process could proceed at tree level 
through bilepton exchange or through various loops. Predictions for the process exist in a variety of 
speculative models including left-right symmetry, 7?-parity- violating supersymmetry, GUT theories and 
several others [pi] . 
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Fig. 17: Proposed beam spill for MECO 



Any possible coupling between muonium and its antiatom will give rise to oscillations between 
them. For atomic s-states with principal quantum number n a splitting of their energy levels 



8Gf G 



MM 



n^-KCLp 



(34) 



is caused, where ao is the Bohr radius of the atom, Gmm is the coupling constant in an effective four- 
fermion interaction and Gf is the weak interaction Fermi coupling constant. For the ground state we have 
5 = 1.5 x 10~ 12 eV x {G M ^/Gp) which corresponds to 519 Hz for G M y^ = G-p. An atomic system 
created at time t = as a pure state of muonium can be expected to be observed in the antimuonium 
state at a later time t with a time dependent probability of 
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where A M = 1/r^ is the muon decay rate (see Fig. |18|). The approximation is valid for a weak coupling 



as suggested by the known experimental limits on G M T7- 




8 9 10 
time [fis] 

Fig. 18: Time dependence of the probability to observe an antimuonium decay for a system which was initially 
in a pure muonium state. The solid line represents the exponential decay of muonium in the absence of a finite 
coupling. The decay probability as antimuonium is given for a coupling strength of G M ^ = IOOOGf by the 
dotted line and for a coupling strength small compared to the muon decay rate (dashed line). In the latter case the 
maximum of the probability is always at about 2 muon lifetimes. Only for strong coupling could several oscillation 
periods be observed. 

The degeneracy of corresponding states in the atom and its anti-atom is removed by external 
magnetic fields which can cause a suppression of the conversion and a reduction of the probability p M y[. 
The influence of an external magnetic field depends on the interaction type of the process. The reduction 
of the conversion probability has been calculated for all possible interaction types as a function of field 
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Fig. 19: The muonium to antimuonium conversion probability depends on external magnetic fields and the 
coupling type. Recent independent calculations were performed by Wong and Hou p4j] and Horikawa and Sasaki 



strength (Fig. 19) [84, 85]. In the case of an observation of the conversion process, the coupling type 
could be revealed by measurements of the conversion probability at two different magnetic field values. 

The conversion process is strongly suppressed for muonium in contact with matter since a transfer 
of the negative muon in antimuonium to any other atom is energetically favoured and breaks up the 
symmetry between muonium and antimuonium by opening up an additional decay channel for the anti- 
atom only [ |86| ] []. Therefore any new sensitive experiment needs to employ muonium atoms in vacuum 

A recent experiment at PSI utilised a powerful signature in which the identification of both con- 
stituents of the anti-tom and their coincident detection was requested after its decay. In this an energetic 
electron appears in the n~ decay. The positron from the atomic shell remains with an average kinetic 
energy of 13.5 eV [Q]. The energetic particle could be observed in a magnetic wire chamber spectrom- 
eter and a position-sensitive microchannel plate (MCP) served as a detector for atomic shell positrons 
onto which these particles could be transported in a magnetic guiding field after post-acceleration in an 
electrostatic device. A clean vertex reconstruction and the observation of annihilation 7's in a pure Csl 
detector surrounding the MCP were required in an event signature [87]. Half a year of actual data-taking 
was carried out at the currently most intense surface muon source, the 7rE5 channel at PSI. The previous 
upper bound on the total conversion probability per muonium atom P M j^ = f p M j^(t)dt was improved 
by more than three orders of magnitude and yielded an upper bound of P M y[ < 8.0 x 10 _11 /5b- Here 
a magnetic field correction Sb is included which accounts for the 0. 1 T magnetic field in the experi- 
mental conditions. It is of order unity and depends on the type of the MM interaction. For an assumed 
effective (V-A) x (V-A)-type four-fermion interaction the quoted result corresponds to an upper limit 
for the coupling constant of G M ^ < 3.0 x lO -3 ^ (90 % C.L.). Several limits on model parameters 
were significantly improved like the mass of the bileptonic gauge boson and some models were strongly 
disfavoured such as a certain Zg model with radiative mass generation and the minimal version of 331 
models [0. 

With a new and intense pulsed beam the characteristic time dependence of the conversion pro- 
cess could be exploited, if only the decay of atoms that have survived several muon lifetimes r M could 
be observed. Whereas all beam-muon-related background decays exponentially, the anti-atom popula- 
tion increases quadratically in time giving the signal an advantage over background, which for a 3-fold 
coincidence signature as in the PSI experiment can be expected to decay with a time constant of r M /3 

4 In gases at atmospheric pressures the conversion probability is approximately five orders of magnitude smaller than in 
vacuum mainly due to scattering of the atoms from gas molecules. In solids the reduction amounts to even 10 orders of 
magnitude. 



29 




T0F - T0F expecte d [HS] 




10 20 



T0F " T0F expected ^ 



Fig. 20: The distribution of the distance of closest approach (Rdca) between a track from an energetic particle in 
the magnetic spectrometer and the back projection of the position on the MCP detector versus the time of flight 
(TOF) of the atomic shell particle for a muonium measurement (left) and for all data recorded within 1290 h of 
data taking while searching for antimuonium (right). One single event falls within 3 standard deviations region 
of the expected TOF and Rdca which is indicated by the ellipse. The events concentrated at early times and low 
i?d ca correspond to a background signal from the allowed decay /j, — ► 2>e2v. In a new experiment such background 
could be suppressed significantly through the characteristically different time evolution of a potential antimuonium 
signal and the background. 



[compare Eq. (35)]. 

Some two orders of magnitude improvement can be envisaged [88] with still no principal back- 
ground arising from the /i — ► 3e2u process or internal Bhabha scattering in which the positron from /i + 
decay would transfer its energy to the electron in the atomic shell and mimic an event which is searched 
for (Fig. |20|). 

The requirements forradiation hardness and rate capability of the set-up are similar to those of a 
\x — > 3e experiment. As before, a common approach to these two measurements may be found. 



3.2.5 Conclusions 

All experiments discussed above would benefit from the highest-possible stop-densities, allowing thin 
targets with minimal distortion of the outgoing particles. In the case of /i + , most promising are subsur- 
face // beams that can be stopped in targets of less than a few mg/cm 2 . 

Better experiments need better detection systems as well. In particular a search for /i — ► e 7 is 
limited by detector technology and it seems unlikely that the sensitivity will rise in proportion to the 
beam intensity. Searches for \i — > 3e and fi + e~ <-> /i~e + may have similar requirements on the set-up. 
High granularity will be required for the tracking system which has to stand the load of the full muon 
decay spectrum. 

From the four tests of lepton flavour conservation discussed here, />e conversion is expected to 
exploit best the increased beam intensity: 

• Intrinsic background can be limited to the decay in orbit which can be kept under control by 
pushing the electron momentum resolution below « 0.5 MeV/c. Care has to be taken to avoid high- 
momentum components in the momentum response function. Ultimately the resolution would be 
given by the uncertainty in the energy loss inside the experimental target. 

• Detector rates can be kept under control thanks to the steeply-falling transverse momentum distri- 
bution. 
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• Prompt beam-correlated background, in particular induced by pions contaminating the beam, can 
be removed using a pulsed proton beam. Alternatively, a cooled muon beam might have a suffi- 
ciently low pion contamination such that only electrons would have to be suppressed. 

In conclusion, all searches would benefit from the projected increase in muon stop rate. In some 
cases the gain in sensitivity over ongoing experiments might, however, be not more than one order of 
magnitude. 

4 FUNDAMENTAL MUON PROPERTIES 

Besides measurements on normal muon decay and the 'classical' searches for the rare decays fi — > erf, 
\i — > eee, \x~ N — > e~N conversion, and muonium-antimuonium conversion (see Section |3~T| ) which are 
strongly motivated by our present thinking towards extensions of the Standard Model, there are a number 
of very interesting experiments on free muons and atomic systems containing muons: 

• precise determinations of the properties characterising the muon (e.g. the mass m«, magnetic 



moment fi^, magnetic anomaly a M , charge and lifetime r M , see Section |2.2[ ), 

• CPT tests from a comparison of and /i + properties, 

• measurements of fundamental constants of nature (e.g. the electromagnetic fine structure constant 
a or the weak interaction Fermi constant Gp), 

• sensitive searches for deviations of muon parameters from Standard Model predictions, including 
novel approaches to search for a permanent muon electric dipole moment edm^, 

• study of nuclear properties in muonic (radioactive) atoms, and 

• applications in condensed matter and life sciences. 

In the past most of these topics have been addressed very successfully and in many cases unique 
information has been extracted. With more intense muon sources they can be accessed with novel and 
much more precise techniques and substantial progress may be achieved. 

4.1 Anomalous magnetic moment of the muon 

The anomalous magnetic moment of the muon, = — 2, has been measured in three experiments 
at CERN, and is presently being measured at BNL [ p0| ] (Fig. 21). A superferric precision magnetic field 



storage ring is employed in which the difference between the spin precession and the cyclotron frequen- 
cies is determined. Currently o M is experimentally known to 1.3 ppm: a^ v = 11 659 202(14)(6) x 10 -10 

& 

The most recent and most accurate theory value in the framework of the Standard Model a M = 
11 659 159.6(6.7) x 10~ 10 (0.57 ppm) appears to differ from the experimental result by Aa M = 43(16) x 
10 -10 which is about 2.6 times the combined experimental and theoretical uncertainties. Although this 
discrepancy may be just caused by an underestimate of the uncertainties in the hadronic contributions to 
the SM value of a^, it is also possible that we are facing hints of physics beyond the Standard Model. At 
this stage, definitive statements cannot yet be made. 

Data already recorded at BNL are expected to provide a statistical error for positive muons which 
is about two times smaller. Data taking with similar accuracy for negative muons is in progress in 
2001, with a potential continuation in 2002. Upon completion, this experiment will either provide more 
confirmed hints of New Physics beyond the Standard Model, or at least limit parameters appearing in 
speculative models. The large number of models to which — 2 has some sensitivity includes super- 
symmetry, lepton compositeness, radiative muon mass models, anomalous W boson parameters, new 
gauge bosons, leptoquark models, exotic fermions and many more [S3]. 
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Fig. 21: A sample of new muon g — 2 data. There is no sign of background after 10 time-dilated muon lifetimes 
of 64 /is. 



As examples of where speculations based on the present difference Sa^ may lead, we note that for 
the constrained minimal supersymmetric extension of the Standard Model (CMSSM) all the preferred 
parameter space is within reach of the LHC, but may not be accessible to the Tevatron collider, nor 
to a first-generation e + e~ linear collider with centre-of-mass energy below 1.2 TeV. In many of these 
models there is a direct connection between a non-standard anomalous muon magnetic moment and other 
processes discussed in this article. For example, in models with muon substructure, a contribution to a M 
from muon compositeness would suggest that the process \i — f could have a branching ratio above 
1CT 13 (see Section [3j]). 

From the experimental point of view, the measurements are predominantly statistics-limited. The 
systematic uncertainties could still be significantly reduced, allowing about one order of magnitude 
higher accuracy before fundamentally new experimental approaches would be needed. At this point, 
the present limitations on magnet and field measurement technology as well as particle pile-up would 
start to impose significant systematic influences on the signal. In order to reduce detector pile-up, a 3.1 
GeV short pulsed (<100 ns) beam of a few xlO Hz repetition rate would be most advantageous for a 
future experiment. A considerable improvement in accuracy could still be envisaged from a new exper- 
imental set-up using a similar approach to the current BNL experiment with a super-ferric storage ring, 
where the main advantage arises from exploiting higher particle fluxes. 

There is, however, a severe limitation to the interpretation of a — 2 measurement, when seeking 
to extract or limit parameters of speculative models, namely that due to the hadronic vacuum polarization. 
Its influence restricts the accuracy of calculations within the Standard Model. Until recently, the hadronic 
contribution was obtained using a dispersion relation and data from electron-positron annihilation into 
hadrons, as well as QCD theory. The accuracy of earlier work [91] was improved significantly by newer 
approaches using hadronic r decay data [89]. It appears that all newer evaluations agree very well with 
the newest and most accurate determination and yield an accuracy of better than 0.6 ppm [92, 93]. Newer 
data for e + e~ annihilation into hadrons are now available from experiments at Novosibirsk and Beijing 
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Fig. 22: Comparison of recent measurements of the muon magnetic anomaly with the most accurate and most 
recent theoretical value JB9||. The difference is 2.6 times the combined theoretical and experimental uncertainty. 

and needs to be included. Further improvements can be expected from the DAPHNE accelerator at 
Frascati, and data on r decays from B factories may also contribute. However, even if the accuracy 
of the theoretical value for hadronic vacuum polarisation could be further improved using even higher- 
quality experimental input, there would still remain an uncertainty arising from the treatment of hadronic 
light-by-light scattering. 

A new and independent g^— 2 experiment would undoubtedly be highly desirable if the experiment- 
theory discrepancy at BNL survives after the experiment is completed. 

4.2 Permanent muon electric dipole moment 

The identification of new sources of CP violation appears to be a crucial requirement for explaining the 
dominance of matter over antimatter in the Universe. Permanent electric dipole moments of fundamental 
particles would violate time reversal (T) invariance and with the assumption of CPT conservation also 
the CP symmetry. 

In the muon g — 2 experiment a hypothetical electric dipole moment of the muon edm^ would 
cause the plane of its spin precession in the storage ring to be tilted against the plane of its orbit. This 
inclination can be measured from an asymmetry of the decay electron distribution with respect to the 
orbital plane. With this method in the latest CERN muon g — 2 experiment a limit for edm^ was 
established at 1.05 x 10~ 18 ecm The ongoing BNL measurements will provide one order of 

magnitude improvement. 



Upper Detector 




Fig. 23: The basic principle of a novel muon edm experiment. /i + 's of velocity v moving in a magnetic storage 
ring with field B find themselves exposed to a motional electric field E oc v x B. In case of a finite EDM the 
muon spin precesses with an angle which grows linearly in time (a), which can be monitored by observing the 
time dependence of the angular distribution of the decay positrons (b). 
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Recently a novel idea was introduced to exploit fully the high motional electric fields in a magnetic 
storage ring. With an additional radial electric field the spin precession due to the magnetic moment 
anomaly can be compensated. In this case an edm^ would express itself as a spin rotation angle out 
of the orbital plane that increases linearly in time (Fig. 23). This can be observed as a time-dependent 
increase of the above/below plane asymmetry. For the possible muon beams at the AGS synchrotron of 
BNL a sensitivity of 10 _24 ecm is expected [94]. In this experiment the possible muon flux is a major 
limitation. For models with nonlinear mass scaling of edm's such an experiment would already be more 
sensitive than present experiments on neutral atoms [95]. An experiment carried out at a more intense 
muon source could provide a significantly more sensitive probe to CP violation in the second generation 
of particles without strangeness. 

It should be noted that the very same experimental technique (and the same experimental set-up) 
could also be utilized for searching for permanent electric dipole moments of a large variety of relativistic 
nuclei [p6|]. Such possibilities were recently explored in a workshop at GSI [ J7|. Results that are superior 
to those extracted from present experiments on neutrons or mercury atoms can be expected. Such a novel 
experiment could be well accommodated at a new and intense ISOL facility associated with the proton 
machine discussed here. 



4.3 CPT tests 

In general, any two measurements of muon parameters for both positive and negative muons constitute 
a CPT test. Such tests were performed for the mass, the magnetic moments and most sensitively for the 
muon magnetic anomaly. In the framework of a rather general ansatz the past muon g-2 experiments at 
CERN have provided the best test of CPT invariance at a level of 2 • 10 -22 which is a more than 3 orders 
of magnitude tighter bound than the mostly quoted K° — K° mass difference and the value extracted from 



measurements of the electron and positron magnetic anomalies Q25|]. From the ongoing measurements 
of a M at BNL one can expect an improvement by at least one order of magnitude in this figure. 

5 BOUND MUON SYSTEMS 
5.1 Muonic atoms 

5.1.1 Muonium 

The hydrogen-like muonium atom (M = /x + e~) is the bound state of a positive muon and an electron 
[pq]. Since both particles may be regarded as structure-less objects with dimensions below 10~ 18 m, 
their electromagnetic interactions can be calculated to much higher precision than for hydrogen isotopes 
containing hadrons. In natural hydrogen, for example, the hyperfine structure has been measured six or- 
ders of magnitude more accurately than can be calculated because of uncertainties in the proton magnetic 
radius and polarizability. For muonium the comparison between theory and experiment of this transition 
could be made with two orders of magnitude higher accuracy than for natural hydrogen. Therefore the 
system gives an important input to the set of adjusted fundamental constants [p9|]. 

Muonium can be produced most efficiently in the atomic Is state. For this reason, precision 
experiments can be earned out on the ground state hyperfine interval and the ls-2s frequency interval. 

Recently an experiment was completed at LAMPF [ ]100| ] in which two Zeeman transitions in the 
ground state were induced by microwave spectroscopy. They both involved a muon spin flip and could be 
detected through a change in the spatial anisotropy of the muon decay positrons. The experiment takes 
advantage of line narrowing that occurs when only such atoms are detected which have coherently inter- 
acted with a microwave field for periods long compared to the average muon lifetime ('old muonium' 



technique [100], Fig. 24). The hyperfine interval was determined to 12 ppb as a test of QED or alter- 



natively as a determination of the fine structure constant a to 58 ppb. The good agreement between this 



a and the value extracted from a measurement of the electron magnetic anomaly (Fig. 25) is generally 



interpreted as the best test of internal consistency of quantum electrodynamics. The comparison between 
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theory and experiment is for the muonium hyperfine structure some two orders of magnitude more pre- 
cise than for electronic hydrogen, where the accuracy of knowledge of the proton magnetic radius and 
polarizability limit the QED calculations more than seven orders of magnitude below the experiment. In 
addition, in the measurements the muon magnetic moment (respectively its mass) were determined to 
120 ppb. 
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Fig. 25: The fine structure constant a has been de- 



termined with various methods (from Ref. [ 101 1) 



Most precise is the value extracted from the mag- 
netic anomaly of the electron. The muonium atom 
offers two different routes which use independent 
sets of fundamental constants. The disagreement 
(the error bars are mostly statistical) seems to indi- 
cate that the uncertainty of h/m e from neutron de 
Broglie wavelength measurements may have been 
underestimated. The degree of agreement between 
the a values from electron g — 2 and muonium hy- 
perfine structure is a measure for the internal con- 
sistency of quantum electrodynamics and the set of 
fundamental constants involved in the evaluation. 



At RAL the ls-2s interval was determined to 4 ppb using Doppler-free two-photon spectroscopy 



with pulsed lasers [102]. The result may be interpreted as a precise muon mass measurement or alterna- 



tively as the by far the best test of the charge equality of muon and electron (2 ppb), i.e. two particles 
from different particle generations. It should be noted that the gauge invariance principle can provide 
charge quantization only within one particle generation. 

The muonium spectroscopy experiments at LAMPF and at RAL were limited largely by statistics 
only. Systematic errors arising mainly from technical imperfections of the pieces of apparatus involved 
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can be expected to be kept under control for one further order of magnitude gain in accuracy for the 
hyperfine structure experiment. The laser experiment can be performed with CW lasers and promises 
some 2-3 orders of magnitude improvement without significant systematics problems, provided 4-5 
orders of magnitude increased muon flux is available. In both cases a pulsed time structure would be 
required. 



5.1.2 Radioactive muonic atoms 

Muonic atoms have been employed both for testing QED vacuum polarisation and to determine nuclear 
parameters, most importantly nuclear mean square charge radii. 

A novel programme is currently being initiated in a collaboration between various universities and 
the Paul Scherrer Institute (PSI) in Switzerland and CERN to pursue the possibility to use muonic X- 
rays to obtain information on very exotic nuclei far from the valley of stability. The exact measurements 
on muonic X-ray spectra can yield the most precise values for the charge radii of nuclei as well as 
other ground-state properties such as moments and even B(E2) transition strengths for even-even nuclei. 
In general, muonic X-rays promise higher accuracy for most of these quantities compared to electron 
scattering. 

The experimental approach proposed is to use ion trapping of exotic isotopes combined with a 
newly developed muon decelerator concept as well as other muon manipulation techniques of PSI. In this 
approach, radioactive ions or isotopes capture muons in their outer atomic orbits, subsequently leading 
to cascade into the inner orbits, resulting in emission of characteristic X-rays of MeV energies due to 
the increase in binding energies for massive muons. Since the capture is of atomic nature, very high 
sensitivities are expected for this approach which, if proven successful, would eventually mean one of 
the greatest new developments in physics of nuclei under extreme conditions. In addition, other physics 
opportunities offered by muon capture in nuclei will be pursued in this collaboration. 

Among the possible scenarios, nested traps for radioactive nuclei and for muons can be envisaged. 



Large formation rates can be expected from a set-up containing a Penning trap Q103[ ] the magnetic field 



of which serves also for a cyclotron muon trap [104]. For muon energies in the range of electron binding 



energies the muon capture cross-sections grow to atomic values, efficient atom production can be ex- 
pected of order 50 systems per second. CERN could be a unique place where such experiments become 
possible. The ideal beam time structure is the shortest possible pulse with non-specific requirements for 
repetition rates. 

Further, nuclear muon capture in both stable and radioactive muonic atoms offers the possibility 
to study weak interactions. Of particular interest is the radiative muon capture in hydrogen, where the 



induced pseudoscalar coupling is a long-standing problem [ 105]. 



5.1.3 Muonic hydrogen 



In muonic hydrogen | |106| ] the Bohr radius of the system is only 100 times larger than the size of the 
proton. Therefore electromagnetic transitions in this system are very sensitive to the size and inner struc- 
ture of the proton. The atom is currently being investigated in an experiment at PSI where the 2s-2p 
Lamb-shift transition will be induced with infrared laser light [107]. This quantity is mostly determined 
by QED effects. A significant contribution arises from the proton mean square charge radius the deter- 
mination of which is the expressed goal of the collaboration at PSI. The proton radius is of significant 
relevance for electronic hydrogen spectroscopy since it gives a large uncertainty in the interpretation of 
the observed ls-2s energy interval. It should be noted that several electron scattering experiments have 
failed to produce a reliable value of the proton mean square charge radius. 

Another important quantity is the hyperfine splitting between the F = and F = 1 states in the 
fip system. A measurement to an accuracy of 10~ 4 or better would yield new information on the proton 
polarisability [108]. Such an experiment can be performed with muonic hydrogen atoms produced by 
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depositing a short pulse of muons in the centre of a small cell in which molecular hydrogen gas at high 
pressure is surrounded by high-Z walls. The average flight time of the system can be measured from 
the timing of muonic X-rays that are emitted upon arrival of the muonic hydrogen atoms at the wall 
where immediate muon transfer occurs. The flight time shortens after laser resonance absorption at 6.8 
/im wavelength due to the immediate collisional quenching and the 0.12 eV gain in kinetic energy from 
recoil in this process. Such an experiment would require an intense pulsed beam with less than 500 ns 
bunch length. If the \i~ beam were polarised the laser transition could be observed as a change in the 



angular distribution of the decay electrons [106]. 



5.2 Condensed matter 

Muon Spin Rotation and Muon Spin Resonance (/^SR) are well-established methods to study bulk con- 
densed matter plQ9p - A typical ^SR experiment requires 10 6 to 10 7 muons. Since any signal is extracted 
from a spin rotation frequency, the resolution of the apparatus is inversely proportional to the duration of 
the muon pulse which therefore should be as short as possible. 

Efforts are currently under way to produce beams of very slow muons which would make it possi- 
ble to expand the [iSR method to research topics such as surface magnetism, surface diffusion, catalysis. 



Bio-science applications are currently being explored Q110Q . They include for example electron transfer 
mechanisms in proteins. Element analysis at surfaces and in thin layers and membranes through muonic 
X-rays can be envisaged. 



For positive muons, moderation to typically 20± 10 eV (Fig. 26) has been demonstrated [111] 
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Fig. 26: The energy distributions for slow muons determined for frozen Ar and N2 targets. 

In frozen noble gas targets the particles have negligible interactions with the moderator atoms once 
they have reached an energy below the materials band gap. The process which generates such a tertiary 
muon beam has some 10 -5 efficiency and requires therefore intense muon sources. The highest flux 
right now is available at PSI [ ]112| ]. 

6 MUON BEAMS 
6.1 Beam requirements 

The experiments discussed in this report require different types of beam. Measurements with stopped 
muons require rather low momentum muons with either a pulsed or a continuous time structure. Other 
experiments would need higher energy pulsed beams (see Table ^). Continuous (or quasi continuous) 
beams would be required particularly for experiments with tight coincidence signatures. 

Low-momentum positive beams can be obtained up to momenta of p s = 29 MeV/c from so-called 



'surface muon' channels [113]. These are beams of muons that originate in the decay of pions that 
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Table 6: Beam requirements for new muon experiments. Given are required sign of charge and the minimum 
of the total usable number of muons J I^dt above which significant progress can be expected in the physical 



interpretation. The experiments which require pulsed beams (see Fig. |27| ) are sensitive to the muon suppression 
Io/Im between pulses of length ST and separation AT. This does not apply (n/a) for continuous beams. Most 
experiments require energies below 4 MeV corresponding to 29 MeV/c momentum. Thin targets and storage ring 
acceptances, demand rather small momentum bites Ap^/p^. 
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T Scenario in which a pulsed beam is utilized. 

* Scenario in which a continuous beam after the muon cooling stage is employed. 
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Fig. 27: A pulsed time structure is important for many stopped muon experiments. Crucial parameters are the 
muon suppression in between the bunches lo/Im, the bunch width ST and separation AT. 



stopped close to the surface of the production target. Their rate increases with p 3 5 up to a maximum 
at p s . Low-momentum negative beams may be obtained from pion decay in the cloud surrounding the 
production target. 

There are important requirements for most experiments on the tolerable momentum bite Ap/p. 
This arises from the fact that either only thin targets are possible in which one would like to stop an as 
high as possible fraction of the beam, or this is determined by the acceptance of subsequent equipment 
like storage rings. Typically one should aim for Ap/p « 1%. Muons outside of this band will not be 
able to contribute to the useful signal, but will most likely contribute to the background. 

Muon beams are regularly contaminated by electrons which originate from -kq decays and pair 
creation. An essential piece of equipment in a low-energy muon beam line is therefore an electromagnetic 
separator (Wien filter). The efficiency of such devices is best for low-momentum spread. (They have also 
an advantage over degrading material, which on the one hand would allow muons to be separated from 
particles of different masses through different fractional energy losses. On the other hand they would 
increase the momentum bite.) 



6.2 Possible realizations of muon beams 

The precise design of muon beams at a neutrino factory complex has not been performed but the fol- 
lowing summarises a number of ways by which one could imagine increasing the muon flux by several 
orders of magnitude over present or foreseen facilities. 

The great strength of the CERN-based Neutrino Factory complex as envisaged in the CERN base- 
line scenario [114], is twofold. Firstly the high intensity offered by the 2.2 GeV/c Superconducting 
Proton Linac (SPL) [115], which is given to provide 4 MW average power in a 75 Hz pulsed mode or 
even possibly 24 MW in continuous wave (CW) mode. Secondly, the great flexibility in the beam timing 
offered by the fact that a proton accumulator and a bunch compressor are foreseen in the complex. These 
rings are necessary to match the basically continuous time structure of the SPL to the pulsed structure 
needed to operate the muon acceleration system as well as the storage ring. 

Figure ^8| shows various locations where targets could be installed along the proton path. 



6.2.1 Continuous beams 

For experiments requiring a continuous beam (// + — > e + 7, fi + — > e + e + e~) two solutions can be en- 
visaged. The benchmark beam for these experiments is the PSI beam for the foreseen [i — > experi- 
ment 



An internal target inside the proton accumulator. This solution, shown in Fig. * s similar to 



that sketched in Ref. [ ]1 16| ] and it would consist in inserting in the proton accumulator a thin target 
(typically 1/4000 of an interaction length), either in place of the H~ stripping foil, or elsewhere in 
the lattice. With respect to the benchmark beam, the gain is obtained from the fact that the beam 
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Fig. 28: Possible locations of muon target along the proton path in the CERN neutrino factory baseline scenario. 



recirculates many times through the target so that essentially all protons - up to tails generated 
by large scatters - end up interacting. This gives right away a factor 20 more intensity than the 
benchmark set-up, where a 5% interaction length target is placed on the proton beam. Since the 
benchmark beam uses a quadrupole capture system, with the typical aperture of 150 mrad, an 
additional gain of a factor 40 or so should be achievable with a more efficient capture system, 
using for instance a solenoidal magnetic field. 

The time structure of the beam is described in Fig. |2^. Because the target is very thin, the protons 
stay in the accumulator for a long time, and the system behaves as an approximately continuous 
muon source. 

This scheme is currently under evaluation. A serious technical problem to be solved is the design 
of the target, which will be subject to very intense heating. This can be in principle solved by 
using a rotating carbon target so that the area of material exposed to the beam (about one square 



centimetre) is removed and radiation-cooled fast enough as in the PSI E and M targets [ 117 ] or in 
a study for the neutrino factory flH8p . Another serious technical problem is the local irradiation 
around the target area. Since all particles produced in the target, except the muons, will have to be 
absorbed locally, this will necessitate a very careful design of the area, which receives the same 
amount of radiation as a beam dump. 

Another possibility would be to use directly the proton beam from the linac operated in high- 
power, CW mode. This solution is less efficient than the previous one, but could constitute a 
back-up possibility. 
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It is therefore seen that improvements due to the target design and capture design, could in princi- 
ple lead to an increase of the flux by almost three orders of magnitude. This will constitute a considerable 
opportunity as well as a tremendous challenge to the design of the experiments themselves. 
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Fig. 29: Sketch of a thin muon production internal target for a quasi-continuous muon beam. 



6.2.2 Pulsed beams 

Muon conversion experiments, fi~N — ► e~N, and various other precision experiments could take sig- 
nificant advantage of a pulsed beam structure (see Table ^ and Fig. ^7|). A pulsed beam allows a time 
delay of many pion lifetimes to reduce the pion contamination. For this rejection to be efficient, a time 
separation of a microsecond between pulses is desirable. Two possibilities come to mind. 

One could choose to place the production target in a proton beam line at the exit of the buncher. 
Taking advantage of the SPL timing flexibility, the bunches could in principle be delivered for this appli- 
cation at a frequency different from the nominal one. 

Another interesting possibility would be to take advantage of the main muon cooling channel, in 
which the particles travel for typically 300 metres, ensuring an efficient clean-up from the pions. The 
muons have a rather broad energy spectrum (typically 200 MeV with an energy spread of 5%) at that 
level, however, and the stopping power may be limited. 

To conclude this section on possible beams, it can be said that the SPL together with its accumu- 
lator and buncher rings offers attractive possibilities. A conceptual design and more precise performance 
estimates of target and beam lines for muon physics at the SPL is now necessary. 

7 CONCLUSIONS 

The main conclusion of this study is that the physics potential of a new slow muon facility, such as the 
one that will become available as a necessary step on the way to building a muon storage ring (Neutrino 
Factory), is very rich and compelling, with a large variety of applications in many fields of basic research. 
Indeed, muon physics, that has already played an important role in establishing the Standard Model of 
particle physics, may provide us with crucial information regarding the theory that lies beyond, proving 
itself to be still far from having exhausted its potential. 
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This new low-energy muon source will have unprecedented intensity, three to four orders of mag- 
nitude larger than presently available. It can have the large degree of flexibility necessary to satisfy the 
requirements of very different experiments, providing muon beams with a wide variety of momenta and 
time structures. Both continuous and pulsed beams are possible. In addition, it is capable of produc- 
ing physics results at the very early stages of the Muon Complex, well before the completion of muon 
cooling, acceleration, and storage sections. 

Only preliminary ideas on the design of this facility are introduced here, suggesting ways by 
which the muon flux could be boosted orders of magnitude above present or foreseen facilities. The 
tasks of detailed conceptual design of target and capture systems and of quantitative estimates of beam 
performances are still entirely ahead of us. It will be now enthusiastically tackled in view of the physics 
rewards that are to be expected. 

A major interest in muon physics lies in the search for rare processes that violate muon number 
conservation. In many extensions of the Standard Model, such as supersymmetry, they may occur at 
rates close to the current experimental bounds. Their discovery would have far-reaching consequences. 
The most interesting processes are n + — ► e +r y, /i + — > e + e~e + , and /x~-e _ conversion in nuclei. We 
emphasise that all the different processes should be pursued. Indeed, the relative rates of the different 
modes provide a powerful tool for discriminating different manifestations of new physics. The muon 
facility discussed here has enough flexibility to allow the study of different muon processes, and promises 
to be more sensitive by at least a few orders of magnitude, when compared with current experiments. 

Experimental prospects on these muon number-violating processes and on the related process of 
muonium-antimuonium conversions can be summarised as follows: 

• The present upper limit on fi + — ► e + 7 (1.2 x 1CT 11 ) is expected to be reduced (or the process 
observed), by 2005, to ~ 10~ 14 . Exploiting the higher intensities available at a neutrino factory 
complex is a nontrivial challenge here, as the main limitation to improving the sensitivity to \x + — > 
e + 7 comes from the severe background of accidental coincidences. A number of paths in detector 
design, aiming at a sensitivity below ~ 10~ 15 , are suggested for further exploration. 

• The reduction in background necessary to improve the sensitivity to \i — ► 3e beyond the current 
upper limit B(fi — > 3e) < 1 x 10~ 12 appears to be an easier task, if detectors working at the high 
beam rates proposed can be developed. Sensitivity to B(/j, —> 3e) > 1 x 10~ 16 or better seems to 
be attainable. 

• The limit on neutrinoless \i~ — > e~ conversion nuclei is expected to be pushed from the present 
6.1 x 10~ 13 down to 1 x 10 -13 soon, while future experiments may be sensitive to conversion 
probabilities down to 1 x 10 -16 . As the experimental sensitivity is completely dominated by the 
performance of the (pulsed) muon beam, this mode should continue to benefit from increased beam 
intensity, and conversion probabilities of 1 x 10~ 18 or smaller could be within reach. 

• Muonium-antimuonium conversion could be further probed by using a pulsed beam and taking 
advantage of the time evolution of the conversion process. One may reach sensitivities to conver- 
sion probabilities as small as 10 -13 , almost three orders of magnitude beyond the present limit of 
8.1 x 10~ n . Detector requirements closely resemble those for fi — ► 3e searches. 

Normal muon decay can also be studied with increased precision. The muon lifetime could be 
measured with a precision which is ten times higher than achieved in the ongoing PSI experiments. 
Improvements on the measurements of decay parameters, via a number of other observables related to 
the muon and electron spin variables, are also possible but are not limited by muon rates, and appear less 
promising in view of the current efforts at TRIUMF and PSI. 

Several other experiments on free muons and atomic systems containing muons also promise sub- 
stantial progress: 
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Table 7: Experiments which could beneficially take advantage of the intense future stopped-muon sources at the 
CERN neutrino factory (NUFACT) 



Type of 


Physics Issues 


Possible 


Previously established 


Present activities 


Projected for 


experiment 




experiments 


accuracy 


(proposed accuracy) 


NUFACT <S CERN 


LldSSlCai 


1 i\i — u — . — ; 

Lepton Number Violation; 


,,~ AT . c — AT 


6 1 X 10 


pet nivin(ii(..t hmi i ~, ^ in — t-V , 

1 ,M. pi OpOseU DlM_ (■) X 1U ) 


in iy 

< 10 




Searches lor Nt'w Physics; 


fi — v erf 


1 2 x 10 — ^ ^ 


piopostU r jl X J-U ) 


< 10 — ^ 


forbidden 


SUSY, L-R Symmetry, 


fi — » eee 


1.0 X 10~ 12 


completed 1985 PSI 


< 10~ 16 


decays 


R-parity violation, 




8.1 X 10 -11 


completed 1999 PSI 


< 10~ 13 


Muon 


G p ; Searches for New Physics; 




IS X 10 _0 


PSI (2x). RAL (1 X 10" 6 ) 


< 10~'' 


decays 


Michel Parameters 


transv. Polariz. 


2 x 10~ 2 


PSI,TRIUMF(5 X 10~ 3 ) 


< 10~ 3 




Standard Model Tests; 










Muon 


New Physics; CPT Tests 


9m - 2 


1.3 X 10~ 6 


BNL(3.5 X 10 -7 ) 


< 10" 7 


moments 


T- resp. CP- Violation 
in 2nd lepton generation 




3.4 X 10~ 19 ecm 


proposed BNL (10 -24 e cm) 


< 5 X 10~ 26 e cm 


Muonium 


Fundamental Constants, f.t,^ L ,a; 


Mhfs 


12 x 10~' J 


completed 1999 LAMPF 


5 X 10~ a 


spectroscopy 


Weak Interactions; Muon Charge 


M ls2s 


1 X 10~ 9 


completed 2000 RAL 


< lo- 11 


Muonic atoms 


Nuclear Charge Radii; 
Weak Interactions 


fi~ atoms 


depends 


PSI, possible CERN 
(< r p >tol0~ 3 ) 


new nuclear 
structure 


Condensed 


surfaces, catalysis 


surface /^SR 


n/a 


PSI, RAL (n/a) 


high rate 


matter 


bio sciences ... 











• A new dedicated muon g — 2 experiment could improve the experimental precision by almost a 
factor 15 compared with present experiments, and by almost a factor of 5 compared with the preci- 
sion expected to be reached by the ongoing BNL experiment (0.35 ppm). This effort may become 
particularly worthwhile if the recently reported discrepancy is confirmed, although hadronic con- 
tributions limit the accuracy of the theoretical prediction. It should also be noted that systematic 
effects on the magnetic field become important at the level of 0. 1 ppm. 

• A search for a permanent muon electric dipole moment (edm) would be a fundamentally novel 
experiment. It could reach a sensitivity similar to or better than present neutron edm searches, 
down to 5 x 10 -26 e cm. 

• Muonium microwave and laser spectroscopy could provide more precise measurements of impor- 
tant fundamental constants like the muon mass, the muon magnetic moment, and the electromag- 
netic fine structure constant. Improvements by one to two orders of magnitude are possible. 

• The combination of intense muon beams and a new ISOL facility offer new physics opportunities. 
Measurements on muonic atoms of radioactive isotopes would reveal properties of nuclei with half- 
lifes down to milliseconds. One can envisage unknown effects in nuclear structure, comparable to 
the recent discovery of halo nuclei. 

• Muonic hydrogen spectroscopy allows measurements of the proton charge radius with much better 
accuracies than what is reached in electron scattering experiments. 

There are also applications in other fields, such as condensed matter physics and life sciences. 
Slow muon beams with four orders of magnitude more flux than the most intense sources at PSI are 
expected to yield new insight in a variety of areas, including surface science, thin films, the chemistry of 
catalysis and dynamical processes in biological molecules. 

We find the discussed experiments theoretically well motivated and capable of providing answers 
to urgent questions in fundamental physics. The described experiments appear technically feasible at the 
projected accuracies (Table ||,f7]). In due time, new international collaborations should be encouraged to 
work on detailed proposals for the possible experiments with forefront detector technology. We enthu- 
siastically recommend to continue efforts towards designing the new stopped muon facility, and we are 
ready to fully support the CERN management in its efforts to achieve such a unique physics opportunity. 
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